Graphene-based photovoltaic cells for near-field thermal energy conversion 
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Thermophotovoltaic devices are energy-conversion systems generating an electric current from the 
thermal photons radiated by a hot body. In far field, the efficiency of these systems is limited by 
the thermodynamic Schockley-Queisser limit corresponding to the case where the source is a black 
body. On the other hand, in near field, the heat flux which can be transferred to a photovoltaic cell 
can be several orders of magnitude larger because of the contribution of evanescent photons. This is 
particularly true when the source supports surface polaritons. Unfortunately, in the infrared where 
these systems operate, the mismatch between the surface-mode frequency and the semiconductor 
gap reduces drastically the potential of this technology. Here we show that graphene-based hybrid 
photovoltaic cells can significantly enhance the generated power paving the way to a promising 
technology for an intensive production of electricity from waste heat. 



A hot body at temperature T radiates an electro- 
magnetic field in its surroundings because of local ther- 
mal fluctuations. In the close vicinity of its surface, at 
distances smaller than the thermal wavelength X t h — 
Hc/(kBT), the electromagnetic energy density is sev- 
eral orders of magnitude larger than in far field [TJ [5]. 
Hence, the near-field thermal radiation associated to non- 
propagating photons which remain confined on the sur- 
face is a potentially important source of energy. By ap- 
proaching a PV cell [5] in proximity of a thermal emit- 
ter, this energy can be extracted by photon tunneling 
toward the cell. Such devices, also called near-field ther- 
mophotovoltaic (NTPV) systems, have been proposed 
ten years ago [1] . In presence of resonant surface modes 
such as surface polaritons, the flux exchanged in near- 
field between source and photodiode drastically over- 
comes the propagative contribution [5H9]. This discov- 
ery has opened new possibilities for the development of 
innovative technologies for nanoscale thermal manage- 
ment, heating-assisted data storage [TO], IR sensing and 
spectroscopy [TTJ [TJ] and has paved the way to a new 
generation of NTPV energy-conversion devices [Tj2 [H] ■ 

Despite its evident interest, several problems still limit 
the technological development of NTPV conversion. The 
main one is the mismatch between the frequency of sur- 
face polaritons supported by the hot source and the gap 
frequency of the cell (typically a semiconductor). Indeed, 
all photons with energy larger than the frequency gap are 
not totally converted into hole-electron pair but a part 
of their energy is dissipated via phonon excitation. Be- 
sides, low-energy photons do not contribute to the pro- 
duction of electricity but are only dissipated into heat 
within the atomic lattice. To overcome this problem, we 
introduce here a relay between the source and the cell to 
make the transport of heat more efficient. Graphene is 
a natural candidate to carry out this function. Indeed, 



this two-dimensional monolayer of carbon atoms which 
has proved to be an extremely surprising material with 
unusual electrical and optical properties |15H17j can be 
tailored by modifying the chemical potential to be res- 
onant between the gap frequency of the semiconductor 
and the resonance frequency of the polariton supported 
by the source. In the context of heat transfer, the role of 
graphene has been recently investigated [HHS] , confirm- 
ing its tunability and paving the way to promising ther- 
mal devices such as thermal transistors. Furthermore, 
a NTPV cell in which a suspended graphene sheet acts 
as source has been recently considered [22]. We propose 
here a modification of the standard NTPV scheme, in 
which the surface of the cell is covered with a graphene 
sheet. As we will show, this enables to exploit at the 
same time the existence of a surface phonon-polariton of 
the source and the tunability of graphene as an efficient 
tool to enhance the source-cell coupling. 
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FIG. 1: Scheme of a typical NTPV cell. A hot source (tem- 
perature T a ) is placed in front of a cell at temperature T c , 
which is typically a p-n junction. The source is heated by 
an external radiation flow, and the temperature of the cell 
is kept constant in time. The radiation flux between source 
and cell is converted into an electrical current inside the cell, 
extracted by means of the two electrodes connected to the 
junction. 
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Figure [T] outlines our novel hybrid graphene- 



2 



semiconductor NTPV system. A hot source made of 
hexagonal Boron Nitride (hBN) at temperature T s — 
450 K, eventually heated by an external primary source, 
is placed in the proximity of a graphene-covered cell made 
of Indium Antimonidc (InSb) at temperature T c = 300 K, 
having gap frequency aj g ~ 2.583 x 10 14 rad s _1 . In the 
frequency domain of interest, the InSb dielectric permit- 
tivity can be described using the simple model [22j 



(1) 




oj > UJ a 



where n r (uj) is the refractive index and the choice ao = 
0.7/im -1 reasonably reproduces the experimental val- 
ues of absorption 23J. As discussed before, it is im- 
portant to choose a source having a strong emission, 
and possibly a surface phonon-polariton resonance, at 
frequencies slightly larger than the gap frequency uj g . 
For our purpose we have chosen hBN which we de- 
scribe (ignoring for the sake of simplicity its anisotropy) 



using a Drude-Lorentz model ei(w) 



>(u> 2 



iTuj)/{uj 2 - wjj + iTu) with = 4.88, lj l = 3.032 x 
10 14 rad s" 1 , uj r = 2.575 x 10 14 rad s" 1 and Y = 1.001 x 
10 12 rad s _1 [13]. This model predicts the existence of a 
surface phonon-polariton resonance at frequency uj spp ~ 
2.960 x 10 14 rad s~ , larger than uj g as wished. 

The hybrid configuration will be compared to the 
case without graphene sheet. The optical properties of 
graphene are accounted for by means of a 2D frequency- 
dependent conductivity (see Appendix for details) [23] 
as already done in the context of a heat-transfer calcula- 
tions in [21] [22] . We are going to study the modifications 
to the radiation flux between the source and the cell, as 
well as to the power produced by the device, due to the 
presence of graphene. In particular, our main scope is to 
investigate whether the presence of graphene is able to 
enhance the efficiency of the NTPV cell and possibly the 
output power as well. 

First, we need a general expression for the heat flux 
exchanged between two planar surfaces. This distance- 
dependent flux tp(d) can be expressed in near-field regime 
as ip(d) = / +00 ^<j>(u,d), where 



4>{oj, d) = Ho: n sc (uj) 



p 



^T p (u>,k,d). (2) 



This expression implies the sum over all the evanescent 
modes of the electromagnetic field (identified by the fre- 
quency oj, the transverse wavevector k = (k x , k y ) and the 
polarization p taking the values p = TE and p = TM) 
of the product of the energy fku carried by each mode 
(w,k,p), the difference n sc (uj) = n{u>,T 8 ) — n(ui,T c ), 
n(ui,T) = [ e hu / k z T - l)- 1 being the distribution func- 
tion inside the reservoir of modes at temperature T, 
and a transmission probability 7^>(w,k, d) through the 
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FIG. 2: Efficiency and ratio of produced electric power in 
presence of graphene. The curves correspond to three differ- 
ent values of chemical potential of graphene. For any distance 
d here represented, the presence of the graphene sheet pro- 
duces an enhancement of efficiency as well as an amplification 
of electric power. The overall effect decreases with respect to 
the source-cell distance d. 



separation gap assuming values between and 1. In 
this Landauer-like decomposition [25-28 the transmis- 
sion probability T p {w, k, d) represents an absolute mea- 
sure of the contribution of a given mode to the energy 
exchange. In the case of two semi-infinite parallel planar 
media this quantity reads }25j 



T p {u,\s.,d) 
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where k z = yjuj 2 /c 2 — k 2 is the z component of the 
wavevector in vacuum, while r\ v and r 2p are the ordi- 
nary vacuum-medium Fresnel coefficients corresponding 
to polarization p and bodies 1 and 2 respectively. For a 
NTPV cell, the expression ([2| has to be modified in order 
to take into account the fact that the cell is a direct-gap 
semiconductor. Hence, the radiative power exchanged 
between the source and the cell is given by 



Piad(d) = / — fcri(u,r s )^(w) 



du) 
2^ 



(4) 



hui n(uj — wq, T c )(j)(ui, d), 



where uiq = eVo/%, Vq being the potential difference at 
which the cell is operating, quantity for which we take 
a value slightly below the theoretical limit Vq = uj g (l — 
T c /T s ) [52]. As for the expression of the electric power 
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which is generated from this flux, it reads 



PMd) 
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fu*j n{u;,T s )4>(uj) 



Hloq n{uj — ujq, T c )(f>(oj, d). 



(5) 



As anticipated, the evolution of both Ppy and the cell 
efficiency i] = Ppy/frad has to be followed when mod- 
ifying a NTPV device, the challenge being their simul- 
taneous enhancement . In Figure [2] we represent both 
the efficiency and the ratio of electric powers for four 
different configurations, namely without graphene and 
with a sheet deposed on the surface of the cell for three 
different values of the chemical potential fi. Regarding 
the efficiency, we see that graphene produces indeed an 
enhancement. For example, at d = 16 nm r\ goes from 
around 10% in absence of graphene to almost 20% for 
fi = 0.5 eV, approaching considerably the ideal Carnot 
limit 1 — T c /T s ~ 33%. We also see that for larger dis- 
tances and for any choice of /i the graphene-modified effi- 
ciencies tend to the one associated to the standard hBN - 
InSb system. As for the electrical power, the presence of 
graphene produces an amplification going up to values of 
the order of 8, showing that both desired conditions are 
met by our modification scheme. 



w spp of hBN. This resonance branch decays with respect 
to the wavevector k 7 and in particular is no longer vis- 
ible around cfc/w spp ~ 100. This cutoff, well-known in 
the theory of radiative heat transfer, is mainly connected 
to the distance between the bodies and is roughly given 
by k, ~ lid [25 1 f26 ] . We remark that the InSb cell 
does not support any surface mode in the frequency re- 
gion of interest. This is no longer true in presence of 
graphene. The modification of the optical properties due 
to the deposed sheet induces the appearance of a reso- 
nant surface mode associated to the graphene-modified 
cell, represented in figure by the dot-dashed line. This is 
not visible for fj, = 0.1 eV, appears for very high wavevec- 
tors for /J, = 0.3 eV, is clearly visible for the largest chosen 
value [i = 0.5 eV. The new physical mechanism generated 
by the presence of graphene can be described in terms of a 
coupling between the two surface modes. For fj, = 0.5 eV 
this manifestly produces an enhanced region of transmis- 
sion probability around the intersection point of the two 
branches, while for /i = 0.1 and 0.3 eV the same coupling, 
taking place at larger values of cfc/w, gives rise to an in- 
crease of the cutoff wavevector k c . The fact that k c is a 
decreasing function of d and that the graphene-induced 
resonant modes exist for high values of the wavevector 
explains why the efficiency enhancement decreases with 
d. 
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FIG. 3: Transmission probability with and without graphene. 
The transmission probability Q is represented in absence of 
graphene, and with graphene for /i = 0.1, 0.3, 0.5 eV. The 
horizontal dashed line represents the frequency oj sp p of the 
surface phonon-polariton resonance of hBN, while the dot- 
dashed lines describe the resonance of the graphene - InSb sys- 
tem. The enhancement of transmission probability is clearly 
visible around the intersection point of the two branches for 
(i — 0.5 eV, whereas the same coupling produces an increase 
of the cutoff wavevector k c for smaller values of the chemical 
potential. 

We now show in Figure [3] the transmission probabil- 
ity 7^(w,k, d) for a source-cell distance of d = 16 nm, 
without and with graphene. We first observe that in ab- 
sence of graphene the modes contributing to the effect 
are concentrated in proximity of the surface resonance 
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FIG. 4: Spectral distribution of radiative flux. The spectral 
flux cj>(u)) is represented for d = 16 nm in absence of graphene, 
and with graphene for /i = 0.1, 0.3, 0.5 eV. The curves show 
an amplification, for any value of the chemical potential, at 
ui = cj S pp. Moreover, for fj, = 0.5 eV the width of the peak 
is significantly enhanced, due to the coupling of surface reso- 
nance modes. 



We finally discuss the effect of graphene on the flux 
spectrum 4>(uj) Q. The first manifest effect is the am- 
plification, for any value of the chemical potential, of 
the peak of the spectrum at u = u; spp . This clearly 
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corresponds to the enhancement of the cutoff wavevec- 
tor k c (i.e. an increased number of participating modes) 
observed for any considered value of /i. Moreover, the 
coupling discussed for fj, = 0.5 eV is here reproduced in 
terms of enlargement of the peak. These curves clearly 
explain that both the radiation exchange and the elec- 
tric power are amplified. Nevertheless, they also allow 
to explain the enhancement of efficiency. As a matter 
of fact, the modifications to the spectral properties are 
more pronounced around w spp , which is larger than tu g . 
Thus, since the region u < u g contributes only to P ra d, 
this amplification results in a higher value of rj as well. 

We have proposed a novel setup for NTPV energy con- 
version in which the PV cell is covered with a graphene 
sheet. The presence of this two-dimensional system mod- 
ifies the optical response of the cell, causing in particular 
the appearance of new surface resonant modes which cou- 
pling to the ones belonging to the source produce a signif- 
icant enhancement of the electric power produced in the 
cell as well as of the overall efficiency of the NTPV cell. 
These graphene-based photovoltaic cells have a consid- 
erable potential to enable ultrahigh-efficiency electricity 
production from thermal energy. 



respectively. We assume the presence of an incoming 
field for z < zq, generating a reflected field in medium 
1 and a transmitted field in medium 2. We then impose 
the continuity of component of the electric field parallel 
to the graphene sheet and connect the discontinuity of 
the magnetic field to the surface current on the sheet, 
proportional to the electric field on the sheet through 
the conductivity a(u>). This procedure gives the following 
expressions of the reflection and transmission coefficients 
for a given couple (w,k) and for the two polarizations 



; (!) ; ( 2 ) I \ 

( 1) ( 2) ' 



tTE 



o i(k^-k™)z _ 



2k 



(i) 



2ik{ 



k^ + k^ + ^oa(cu)uj 
f (,Ak {1) f f(1 u(2) , 



(1)1.(2) ' 



ea ( w )*«+ ei ( W )*W + 2^ 



J(k^-k^)z _ 



2ei(w)fci 2) 



e 2 (u J )k^ + e^ + ^^ 2} 



(7) 



Appendix 



where 



We give here a brief overview of the optical properties 
of graphene used throughout the calculation presented 
in the paper. The response of the graphene sheet is de- 
scribed in terms of a 2D frequency-dependent conductiv- 
ity cr(w), written as a sum of an intraband (Drude) and 
an interband contribution respectively given by |24j 
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where G{x) = sxah(x/k B T)/[cosh(fj,/k B T) + 
cosh(x/k B T)]. The conductivity also depends on 
the temperature T of the graphene sheet (assumed equal 
to the temperature of the cell) , on the chemical potential 
(i and on the relaxation time t, for which we have chosen 
the value r = fCT 13 s [29]. 

In order to deduce the heat flux between the source 
and the graphene-covered cell we need the expression of 
the reflection coefficients for an arbitrary frequency uj, 
wavevector k and polarization p. To this aim we assume 
that the graphene sheet is located in z = z and that 
it separates two non-magnetic media 1 (z < zq) and 2 
(z > zq) having dielectric permittivities £i(w) and s 2 (oj) 



kf = v^H 



k 2 
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is the z component of the wavevector inside medium i 
(i = 1,2). For the purpose of our calculation we simply 
have to take £\{ui) = 1, obtaining 
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